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Abstract 

In this study we show that the sum of the powers of arbitrary products 
of quantum spin operators such as {S^)\S^)™'{S^)^ can be reduced by one 
unit, if this sum is equal to 25" + 1, S being the spin quantum number. We 
emphasize that by a repeated apphcation of this procedure all arbitrary spin 
operator products with a sum of powers larger than 25 can be replaced by 
a combination of spin operators with a maximum sum of powers not larger 
than 23. This transformation is exact. All spin operators must belong to 
the same lattice site. By use of this procedure the consideration of single- 
ion anisotropics and the investigation of the magnetic reorientation within a 
Green's function theory are facilitated. Furthermore, it may be useful for the 
study of time dependent magnetic properties within the ultrashort (fsec) time 
domain. 

PACS: 05.50.+q, 75.10.Jm, 75.30.Gw 



A renewed interest in the properties of spin operators and their products is caused by 
the investigation of the magnetic anisotropy and other magnetic properties of thin films and 
multilayers [^. In general the lattice anisotropics are defined by the expansion of the free 
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energy F{T, Q, $) in powers of the direction cosines, or by a complete set of spherical har- 
monics Yim{Q, 0, with O and $ being the polar and the azimuthal angles, respectively. 
The coeffi cents of this expansion are temperature dependent, these effective anisotropy co- 
efficients are measurable quantities To compare the measurements with first principle 
electronic calculations of the microscopic anisotropy parameters at T = |^,|^ , the temper- 
ature dependences of the effective anisotropics have to be known [0-^. 

The magnetic lattice anisotropy is mainly a single-ion property. Corresponding to the 
atomic lattice symmetry, they can be represented in the Hamiltonian by powers of the 
components of the three-dimensional Heisenberg spin S"' , a = x,y, z, and mixed products of 
these spin operators. By use of this generalized Heisenberg-type Hamiltonian the magnetic 
properties, in particular the effective anisotropics, can be determined by the usual rules of 
Statistical Mechanics |0,p,|lO|,|ll| . Note that in addition to the diagonal operators such as 
(S*^)" or (S^S^)"', also non- diagonal operators may occur. Here we have introduced the non- 
hermitean raising and lowering spin operators 5*^ = S^±i S"^, assuming the quantization axis 
along the z- direction. This direction may be the easy axis of the magnetization determined 
by the anisotropics, for example. Due to an external magnetic field the spins are tilted away 
from the easy axis. Then by application of the usual rotation properties of spin operators 



1^ , non-diagonal spin operators are introduced in the Hamiltonian. 

For a single-spin problem, the respective expectation values of spin operators can be 



easily deduced from simple diagonalization p,10,ll|. However, in the case of interacting 



extended systems, in which the magnetic moments are coupled by the exchange interac- 
tion, approximations are necessary. The exchange coupling is in general much stronger than 
the involved anisotropics. Within a mean field approximation the system is reduced to an 
effective single-particle system, subject to a molecular field due to the exchange coupling 



P JTO|JTl[] . This approach yields a bad discription of the magnetization in particular for low- 
dimensional systems, since collective excitations (spin-waves) of the Heisenberg magnet are 
neglected. These spin fiuctuations are considered by spin-wave theories. For instance, the 
Holstein-Primakoff approach has been applied frequently 0,0], which is valid for low tem- 



peratures. On the other hand, the Green's function method |T5|-[T7[] and the Schwinger-Boson 



approximations |T8| , p!9| can be used also at higher temperatures, since they consider partly 
interactions between spin-waves. The Green's function method is a hierachical approach, 
i.e. in the respective equation of motion appear higher order Green's functions consisting of 



a larger number of spin operators |TB,[TB|. In order to obtain a closed set of equations for the 



Green's functions, the higher order ones have to be decoupled, i.e. approximated by lower 
order ones at a certain level of hierarchy. 

The single-ion lattice anisotropics are taken into account by local Green's functions, i.e. 
consisting only spin operators of a single lattice site. Such local Green's functions with 
non-diagonal operators appear within the spin wave treatment of the magnetic reorientation 



r?| , or in the case of a noncollinear magnetism. In this study we will show that under 
certain conditions the higher order local Green's functions appearing in the equation of 
motion, or the respective expectation values as obtained from the spectral theorem, can be 
replaced by lower order ones. This will be the case for general products of spin operators 
such as or if the sum of their powers (the 'total power') 

is / + m + n = 2S + 1, S being the (integer or half- integer) spin quantum number. In this case 
the total power can be reduced by one unit, resulting then in lower order Green's functions 
or expectation values. Any combination of spin operators can be reduced, in particular also 
non-diagonal ones such as {S^)"^{S^)'^^~^^~"^ . This procedure is exact. The reason for this 
property is the fixed magnitude of the spin operator S, i.e. (S)^ = {S^Y + (S^)^ + (S^)^ = 
S{S + 1), thus the three components are not independent of each other. Furthermore, by 
a repeated application of this procedure also spin operator products with a sum of powers 
larger than 25" + 1 can be reduced correspondingly to a total power equal to 2S. We are not 
aware in the literature about such a reduction in general. The outlined procedure seems to 
be a special case of the relationships between irreducible unit tensor operators with different 
rank p!2|j20| . In a forthcoming paper we will investigate this connection into greater detail. 

Quantum mechanical spins are characterized by the commutation relations [S'+,S'^] = 
2S^ and [S""^, S"^] = ±5"^, in units of h = 1. As mentioned, the spin operators under 



consideration belong to the same Heisenberg spin S, thus a particular lattice site index is 
omitted here. In case of time dependent spin operators the same time for all operators is 
assumed. The description of the reduction of the spin operator products is performed within 
several steps. As an example we show the respective procedure of the particular product 
{S^)"^{S^y^^^^"^. Other combinations can be obtained similarly, or by application of the 
commutation relations between spin operators. 

1) First we state the general property for the raising and lowering operators, 

(5-)2^+l = (5+)2^+l = 0. (1) 



2) From the the general identity ||T3|,[TB| 



fliS'-^) = 0, (2) 

i=-S 

one can deduce the first reduction in the total power of spin operator products, namely 

2S 
i=0 



{sr'^' = j2o^f\sr, (3) 



(S) 

the ai being rational coefficients dependent on 5", which can be determined for each 5* 
with the help of a recursion relation from coefficients of lower S. In particular one obtains 

5-1 

integer S : {S^''^' = E «2? (ST , 

S-l/2 

half - integer S : (5^)=^^+^ = ^ "2? (Si'"' ■ (4) 

i=0 

(S) 

The coefficients for the minimal and the maximal index i can be given in closed form, 

5-1 

integer S : = -S^ 4^"'^ = (-1)"^-' H - , 



i=0 

S 

(5-1) 



5-1/2 

half - integer S : a^, = -S^ a^^"'^ = {-if-'/^ U - 



i=0 

5-1/2 

(5-1) 



'^25-1 — '^25-3 

i=0 



For the other indices the recursion relations are given by 

cyi = cyi~^\-S'c,ir\ (6) 

with 2i = 2, . . . ,2S — 4 for integer S, and 2i = 2, . . . ,2S — 3 for half-integer S. 

3) Consider now the general spin identity 

[S'Y = S{S + 1) - S' - S-S+ . (7) 

For the determination of the products (^5*-)"^ (^5*^)25+1-771 j^iultiply Eq.(0) from the right 
subsequently by (5*^)', / = 2S' — 1, . . . , 1. Calculate at first for / = 25' — 1 with the help of 
Eq.(^ the product 

2S-1 

(S'fs ^ J2 '2^-1) {SJ - S-S+ {S'f^~^ , (8) 

1=0 

the coefficients f]^^''^^~^^ depend on the a\^^ and have to be determined for each S. Calculate 
then subsequently for / = 25' — 2, 25* — 3, . . . , 1 the products 

I 25-1 

(sT' = E f^f''\sy + S-S+ E , (9) 

i=0 i=l 

(s n (s n (s 2,S—1) 

with corresponding coefficients f3l ' and 7^ ' dependent on 5', with 72s'_i = — 1- Within 

this procedure the product (5*^)'+^ appearing for each value of / has to be replaced by the 

result of the preceeding calculation for Z + 1. Finally, for / = one obtains an expression for 

5"^ given by 

25-1 

= + 5-5+ E ■ (10) 

Eqs-d^), (|), (y), and (|T0|) form a set of equations for the powers of 5'^ ranging from 
to 25 + 1. The products 5^5+(5^)* remain unchanged at present. These equations can 
be viewed to represent relations between spin operators as well as between the respective 
expectation values. 

4) By application of Eqs. - (JT^) we are now in the situation to reduce the sum of powers 
of the spin operator products (5^)^(5'^)^'^^^"™' by one unit. Consider at first the product 
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{S-y^S' (m = 2S), and replace the operator S' by Eq.(|TO|). Due to {S-)^^+^ = 0, Eq.(|), 
the terms consisting S^S^ on the right hand side of Eq.(p!0|) vanish, yielding 

{s-y^s^ = 5f {s-f' = °) {s-y' = s {s-f' . (11) 

The latter identity is obtained from inspection. The total power of the product {S~f^S^ 
has thus reduced by one unit. 

5) Consider now subsequently the products (^5'-)™-(5'z^2S'+i-m ^j^]-^ m = 2S — 1, . . . ,1. As in 
the preceeding step, replace the operator (^5'2)25+i-m, ^\^q j-^gjp Qf Eq.(|^) for / + 1 = 25 + 
1 — m. Use then the result for the preceeding spin operator product (S-)"^' {S^f^+^'"'' with 
m' = m+1. Within the course of these calculations one has to commute the product S^{S^y 
by use of S+{S'y = {S' - iyS+, replace the product S'S+ by the identity Eq.(0), and 
finally solve the result for {S^)"^ {S^Y^'^^~"^ . Together with Eqs.@,(|llD a set of equations 
is obtained for the products 

2S-m 
1=0 

for < m < 2S. As can be seen, the products of the left hand side of Eq. ([1^) , with the 
total power equal to 25* + 1, have been replaced by a sum of spin operator products whose 
maximal sum of powers is 25". Thus the total powers of the original spin operator products 
have been reduced by one unit. The rational coefficients S^^'"^^ depend on the spin quantum 
number S. For m = they are given by the coefficients ap'^^ of Eq.(^. In Table 1 we 
present results for ^j"^'™-* for several values of the spin quantum number S. For the maximal 
index i = 2S — m one obtains from inspection 525!!^ = m{2S + 1 — m)/2, furthermore for 
integer S one has 5g = 0. As an example for the reduction procedure and for the use of 
the Table 1 we show exphcitely the respective results for the case 5 = 2. The products of 
spin operators (S'~)'"(S'^)^'^^^^™ with a total power of 25* + 1 = 5 is replaced by a series 
with increasing powers of S"^, and a maximum total power of 2S = 4. 

(5^)5 = _4 S' + 5{S'f 

S~ [S^y = -2 + S' {S'f + 2 S-{S')^ 
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{S-y {S')^ = -2 {S-f S' + 3 {S-f {S'f 
(S- f (S' f = -2 (5-)^ + 3 (S-f 

{S-f = 

6) To consider the spin operator products {S^)'^''^~^^^"^{S^)"^ , Eq.(|^) has to be multiphed from 
the left by {S'Y. Eqs.® - (|lO|) remain unchanged if S'S+iS')' is substituted by {S^YS-S^, 
since S^S^ is diagonal in the representation. Thus the products (5'^^25+i-mj^^+-jm ^^^^ 
replaced by 

2S-m 

(^S'fS+l-m^g+^m ^ ^ ^g+^m _ (^g) 

i=0 

7) The corresponding reduction of spin operator products with an arbitrary order of the 
operators S^, S~ , and S^, for example (S^)'^^^^^'^ (S^)"^ , can be deduced from Eqs.(|12]),(|13D 
by the usual commutation relations of the spin operators. Products consisting both 5*+ and 

may be facilitated with the help of Eq.(|^). If the total power is larger than 25 + 1, at 
first a particular subset of spin operators can be reduced. Then by repeated application of 
the procedure the complete spin operator product is replaced by a sum of products with a 
maximal total power equal to 25*. 

8) At the end of this derivation we like to remark a few important points. 

(i) The described procedure for the reduction of the local spin operator products is exact, 
no approximation has been used. Products of spin operators located on different lattice 
sites, for example due to the exchange interaction, cannot be considered by the procedure 
outlined here. 

(ii) The above results are valid for a sum of powers equal to or larger than 25+1. The 
operators (S*^)™ in Eqs.(|12]),(|T3|) must not be cancelled, since this will lead to wrong results. 
For instance, cancellation of the common factor {S~)"^ in Eq.(|T2|) leads to the relation 
between powers of given by (5*2^25+1-™ _ 5^'^'™-' [S^y. However, this relation is 
only correct for m = 0, i.e. Eq.(RD ! 



(iii) Facilitating the powers {S^Y^^ by replacing the products S~S~^ on the right hand side 
of Eq.(^) - (|10|) with the help of Eq.(|^) does not help since this leads to identities. 

(iv) Similar as for the coefficients a in Eq.(^ there should exist also recursion relations for 
the coefficients /3, 7, and 6 of Eqs.(H) - (|13|), connecting respective coefficients for different 
spin quantum numbers S. We have not taken here the effort to deduce them. 

In conclusion, in this study we have shown how the total power of arbitrary products 
of spin operators {S'^Y {S~)"^ (S^)"' can be reduced by one unit, if the sum of its powers is 
equal to 25* + 1. In addition we emphasize that by a repeated application of the procedure 
all arbitrary spin operator products with a total power larger than 25 + 1 can be replaced 
by a sum of spin operators with a maximum total power of 2S. These transformations 
are exact. It is important to note that the procedure outlined here is valid for local spin 
operator products only, i.e. all involved spin operators belong to the same lattice site. For 
example, the magnetic anisotropics, determining the direction of the magnetization, are 
usually single-ion quantities. We note that the consideration of such a single-ion anisotropy 
of l-th order are only supported by spin quantum numbers 25* > / |]TD|,[T2|,pD|,pT]| . Vice 



versa, for a given S the maximum order of distinct anisotropy terms as considered in the 
Heisenberg Hamiltonian is limited since anisotropy orders larger than 2S can be replaced 
by lower order ones. This is due to exactly the same circumstances leading to the reduction 
of spin operator products outlined here. Note also that single-ion anisotropy terms can 



be replaced by irreducible tensor operators [10,12,PQ,piI|. Furthermore, the investigation of 



the magnetic properties within a Green's function theory can be facilitated |T4|-[17[|. For 
instance, the higher order Green's functions appearing in the equation of motion can be 
replaced by lower order ones without approximation, if the involved local spin operators are 
reduced along the lines described here. 

Localized, quantum mechanical spins are present in the strongly correlated rare earth 
magnets, in which the localized 4/- electrons cause the magnetic state. On the other hand, 
the 3d- magnets are itinerant electron systems, the electrons causing the magnetism are not 
as localized as in the 4/- systems, and the quantum mechanical nature of the equilibrium 
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state of the magnetic properties may be hidden. This can be seen for example from the 
arbitrary values of the resulting magnetic moments. Consequently, its magnetic properties 
are usually represented by three-dimensional classical spin vectors S —>■ oo. However, the 
quantum mechanical nature of the magnetism of these band magnets should be visible in the 



ultrashort time domain, for example within the magnetic relaxation after an excitation E 



Such ultrashort temporal properties in the fsec- time domain can nowadays be investigated 
experimentally ||23|j2^ . For the theoretical study of these time dependent magnetic prop- 
erties, which can be accomplished e.g. by application of a time dependent Green's function 
(Keldysh-) formalism ||25|, the reduction of the spin operator products may be a very useful 
tool. 
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s 



m 



(S,0) j^{S,0) j^iSfi) ;^{S,0) ;^{S,0) ^(S,0) 



ir{0,U) ir(0,Vj r^o.u; rVO.Uj r 
Oq Oi O2 O3 0, 



.(5,1) .(5,1) .(5,1) .(5,1) .(5,1) 
Oq O2 O3 O4 



(5,2) ;^{5,2) ;^(5,2) ^(5,2) 



Oq 0]_ 0. 



m 



^(5,3) ^(5,3) ^(5,3) 



(5,4) .(5,4) 



Oq Oi 



9_ 
'16 



225 
64 





1 


-4 




10 

4 



259 
■ 16 





5 




35 
4 



45 

32 









1 






1 

4 


3 
2 




-2 


1 


2 


9 

16 


25 
4 


5 
2 



15 
16 



2 
-2 
-1 



3 
2 

-2 

15 



23 
■ 4 



15 
■ 4 



Table 1: Coefficients (5- '"^ occurring in Eqs.(12),(13) for several values of the spin quantum number S = 1/2, 1,3/2,2,5/2. 
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